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I. INTRODUCTION 



In recent years there has been significant progress in the study of the spectrum of hadrons, 
as well as their non-perturbative structure and form factors, through approaches that are 
manifestly covariant and which accommodate both dynamical chiral symmetry breaking 
(DCSB) and quark confinement Covariance provides efficient and unambiguous access 

h h fl 

to form factors [2, l3|, |J] and parton distribution functions [5]. The covariant approach 
to modeling QCD with which we will be concerned is based upon the Dyson-Schwinger 
equations (DSE) of QCD [1] . These are coupled integral equations of motion for the 
QCD Green's functions or n-point functions. The DSE for any particular n-point function 
involves higher (n+l)-point functions and so on. One needs to employ a truncation of the 
coupled equations for practical calculations. 

One of the simplest and most studied DSEs is the quark propagator DSE, the quark gap 
equation 

S-\p) = Z, S^\p) + C^Z, g'D^^ip - g) 7^ X S{q)r,{q,p), (1) 

Jq 

where Sq^{p) = i'j ■ p + rribm, fntm is the bare current quark mass, T^(q,p) is the dressed 
quark-gluon vertex function, = d'^q/ {2ttY denotes a loop integral regularized in a 
translationally-invariant manner at mass-scale A. Here A^) is the vertex renormaliza- 

tion constant to ensure To- = 70- at the renormalization scale /i and Z2(/i^,A^) is the quark 
wave function renormalization constant. The color factor in Eq. ([1]) is Cp = {N^ — l)/2Nc- 
The general form for S{p)^^ is 

S{p)-'=^^■pA{p\^^'') + B{p\^^') (2) 

and the renormalization condition at scale p^ = is S{p)^^ —>■ i'j ■ p + Tn{fi), where m(/i) 
is the renormalized current quark mass. One can see from Eq. ([1]) that, in principle, one 
must first solve the DSE for the gluon 2-point function as well as the DSE for the dressed 
quark-gluon vertex. In the renormalization group improved rainbow approximation Q], one 
replaces the combination of those two objects by bare vertex times an effective gluon 2-point 
function. The latter has a phenomenological infrared part that joins smoothly to the 1-loop 
renormalization group result for the ultraviolet running coupling. Typically one infrared 
parameter is sufficient to generate the empirical value of the quark condensate. When the 



same kernel and dressed propagators are used for the Bethe-Salpeter equation (BSE), this 
ladder-rainbow (LR) model gives a very good account of light quark vector and pseudoscalar 
ground state mesons [7] and many other observables - see Ref. [l| for a review. 

As this LR truncation is a non-perturbative truncation, one does not have an easy handle 
on its accuracy. Recent efforts have been made to perform calculations of meson masses 
beyond ladder-rainbow truncation 

aai 

|lO|. These studies employ a quark-gluon vertex 



that is dressed by certain classes of diagrams constructed from 2-point gluon functions 
alone. The results of a truncation back to a bare vertex indicate that the error in the LR 
truncation is typically 10-30 percent, depending on the class of vertex diagram that is taken 
into account. To make such calculations feasible, the gluon propagator is simplified to a 
momentum delta function so that non-linear integral equations are reduced to non-linear 
algebraic equations. The sum of ladder-rainbow diagrams for the vertex was implemented 



in Ref. [8|, while in Ref. [9| the important attractive effects of dressing with 3-point gluon 
functions were added phenomenologically. The 4-point gluon functions have so far been 
neglected altogether. In our previous work [l^, we studied an improved ladder-rainbow 
summation in which internal vertices were self-consistently dressed. In the current work we 
explore a scheme for including all possible diagrams involving vertex dressing with 2-point 
gluon functions and explore the effects of this generalization on the behavior of the quark 
propagator amplitudes. This entails the inclusion of crossed gluon line diagrams. A study 
of fully dressed quark-gluon vertex in scalar ^p^x theory with Munczek-Nemirowski gluon 
propagator of Ref. showed that the solutions for quark propagator functions show large 
variations in the region of the small momenta in comparison with those corresponding to 
LR truncation of the quark DSE. 

To demonstrate the relevance of the current dressing scheme, consider the vertex with 
exactly 2 gluon lines. Fig. [1] In our previous work with the most general vertex dressing 
scheme to date [Iff] , we extended the ladder-summed approximation of Refs. 0, (Fig. [TK) 
to self-consistently include the dressing of internal vertices (Figs. Wp and [lb). The same 
procedure was extended to a general number of gluon lines and the contributions were 
summed until convergence was reached for this improved ladder-rainbow sum. Crossed 
gluon line diagrams, such as Fig. [T]i were not included. Convergence in this improved 
ladder-rainbow sum was found in the vertex itself, in the quark DSE solutions, and in the 
pseudoscalar and vector meson masses. 
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The diagrams in the improved ladder-rainbow sum have the same simple color structure 
and the associated color factors are powers of one basic color factor. The color structure 
analysis of the diagrams omitted from the improved ladder-rainbow class shows that they 
can have significantly larger color-factors and that, at the very least, their role should be 
explored. For example, with A^^^ colors, the color factor associated with Figs. [T^, [Dd and [T]3 
are 1/{2N^'^, whereas for the crossed gluon diagram of Fig. [T]i the color factor is {N"^ + 
l)/(2Ai"c)^. In general, similar examples occur at higher order. Also working to magnify 
the contributions to the vertex from many gluon lines is the fact that the number of such 
diagrams increases rapidly. A formal analysis of these additional diagrams in terms of l/A'c 



expansion is presented in Ref. 



121. 




FIG. 1: (Color Online) Quark-gluon vertex diagrams at O(g^), neglecting 3-gluon and 4-gluon 
coupling. The large red circle denotes the second order vertex and the small black circles indicate 

WW 

dressed propagators. Diagram (a) is part of the ladder sum [a, |9(, while diagrams (b) and (c) are 



added in the improved ladder sum 



10|, and diagram (d) enters in the extension considered in this 



work. Crossed gluon diagrams, such as (d), prove to be important because of their significantly 
larger color factors. 



II. THE FULLY DRESSED VERTEX WITH TWO-POINT GLUON LINES 

We employ Landau gauge and a Euclidean metric, with: {7^,71/} = 25^y] 7^^ = 7^; 
and a ■ b = Ylt=i '^^^ dressed quark-gluon vertex for gluon momentum k and quark 
momentum p can be written igt^Tcj{p + k^p), where = and A'^ is an SU(3) color 
matrix. In general, r„{p+k,p) has 12 independent invariant amplitudes. We are particularly 
concerned in this work with the vertex at = 0, in which case the general form is 

r^(p) = ai(p^)7^ + a2(/)7 ■ P Pa - a:i{p^)ip„ + ai{p^)i'-f„ j ■ p (3) 
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where are invariant amplitudes. 

We propose a scheme for dressing the vertex function that includes all possible contri- 
butions constructed with two-point gluon lines up to n gluon lines. We can decompose the 
vertex function into a sum of all contributions with exactly i gluon lines: 

= E (4) 

where r° = Zi7^. 

In previous work we developed a scheme to include all vertex diagrams that have a ladder 



structure but further improved it by the self-consistent dressing of all internal vertices |10 |. 
In that case the contribution is generated from three vertex contributions having a smaller 
number of gluon lines by adding one gluon ladder rung with dressed vertices. If the number 
of gluon lines in the three vertex contributions are denoted j, k and /, then summation is 
made over j, k and / such that j + A; + / + 1 = i. That iterative self-consistent ladder scheme 
is described by 



i=j+k+l+l 

X Tiip +k,q + k)Siq + k)Tliq + k, g)S(g)r^(g, p), (5) 

for i > 1. Here, and in all vertex constructions that we discuss and implement, the quark 
propagators that are internal to the vertex are self-consistently dressed using the gap equa- 
tion, Eq. ([1]), with the dressed vertex at the order being discussed. The color factors of all 
contributions included in this improved ladder vertex are correctly generated by the iter- 
ations of implied by the above. Such a vertex, generated by 2-point gluon lines, is 
repulsive and as a consequence one finds that the dominant vertex amplitude < 1. 

However, the available lattice data indicate that Q;i(p^) > 1 much as it is in an Abelian 
theory. Certainly in the ultraviolet, where a one-loop analysis is reliable [l^, ai(p^) — ^ 1 
from above because of the strong attraction is provided by the 3-gluon coupling term with 
its dominant and large color factor In Refs. jo], [l^, it was argued that an effective way 
to simulate such physics in studies like the present is to assume the momentum dependence 
of the corresponding 3-gluon coupling contribution to the vertex is similar enough to Eq. ([5]) 
to allow the color factor there to be replaced by the sum of the two, which is Cp, the color 
factor of the quark DSE in Eq. ([1]). To reproduce lattice vertex data, and to explore the effect 




FIG. 2: (Color Online) Quark self-energy type contribution in dressing the quark-gluon vertex. 



of 3-gluon coupling in a slightly more general way, Refs. [y, [lOj, employed Eq. ([5]) with color 



which corresponds to omitting all consideration of 3-gluon coupling. Subsequently, with the 
choice C = 0.5, we explore the likely influence of 3-gluon coupling with the present vertex. 

In this present work we explore the addition of the remaining class of vertex diagrams 
that are beyond those in the improved ladder sum outlined above, but still obtainable in 
terms of 2-point gluon functions. 

To generate such a vertex, we use a recursive algorithm for obtaining all the contributions 
with exactly i gluon lines, , from the diagrams contributing to F^^^. We take each diagram 
contributing to F^~^ and consider the diagrams produced by making all possible insertions 
of a single gluon line so that it starts and ends on a fermion line without being simply a 
fermion self-energy insertion (see Fig. ([2])). Then, to avoid double counting, a check is made 
as to whether the resulting diagram was already produced from any previous construction 
at the same order. One can easily see that the algorithm works for constructing the one-loop 
diagram from the bare vertex and from that the two- loop diagrams of Fig. ([T]). The steps 
of the algorithm are illustrated in Fig. ([3]). 

It is trivial to prove, using mathematical induction, that this algorithm produces all 
possible diagrams involving n gluon lines in dressing the vertex for any n > 1. In the previous 
paragraph we demonstrated that all the first and second order diagrams are generated using 
this algorithm. Let us suppose, for some n > 1, that we start with the correct diagrammatic 
content of but that the algorithm for failed to produce a certain valid diagram 
with n gluon lines. It is easy to see that we can always find one gluon line, which can be 




factor CCp where C varies in the range 



1 



to 1 . Here we will initially use the lower limit 



2Nc Cp 
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FIG. 3: (Color Online) Generating all the second order diagrams from the only first order diagram 
by making all possible insertions of a single gluon line into the first order vertex diagram. The 
letters indicate the starting point of the inserted line. We note that there are several redundant 
diagrams generated, which must be removed. 

removed from that diagram so that the remaining n — 1 gluon lines are not simply quark 
self-energy insertions. Hence the produced parent vertex is a valid member of and the 
supposition that all possible insertions into it are not capable of generating all members of 
is incorrect. 

In fact, the gluon line that can be removed is easy to find in a diagram. If a gluon line 
attached to the end of a quark line cannot be so removed, then the other end of this gluon 
line separates ends of a second gluon line that would otherwise be a self-energy insertion. 
This second gluon line is the one that can be removed. 



7 




FIG. 4: (Color Online) Constructing the set of numbers uniquely identifying one of the second 
order diagrams. 



n n 

( 1^2, 1, 0, n-^ ) 

FIG. 5: (Color Online) The sets representing all diagrams with exactly n gluon lines are obtained 
from those corresponding to diagrams with exactly n — 1 lines by inserting a pair of n's into the 
set and making all possible permutations, so they are never next to each other. 



The algorithm described above is easy to implement in practice by constructing a set 
of numbers uniquely identifying each diagram. We build the set by enumerating the bare 
quark-gluon vertices in a diagram with n — 1 gluon lines from 1 to — 1 and assigning the 
same numbers to the vertices attached to the same gluon propagators. We assign to the 
external gluon vertex. An example of such construction is shown in Fig. (jl]). 

To construct the vertices with n gluon lines we insert a pair of integers n into the set 
described above, so that they are not next to each other, as illustrated in Fig. ([5]). We relabel 
the resulting set in ascending order and check whether the final set was already generated. 



1. Algebraic Analysis 

In order to solve the gap equation with our dressed quark-gluon vertex one needs the gluon 
2-point function. As in the previous studies of this kind, we again employ the Munczek- 
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Nemirovsky Ansatz [14] for the interaction kernel in the form 

g'D,^{k) - ^) (2^)'^'5'(fc), (6) 

where the parameter is a measure of the integrated kernel strength. This simplifica- 
tion reduces the mult i- dimensional integral equations for the vertex and the dressed quark 
propagator to algebraic ones that make the analysis feasible. 

The algebraic form of the gap equation for the quark propagator corresponding to the 
gluon propagator in Eq. ([6]) is 

S"^{p) = i-i ■ p + rribrn + G'^l^^S{p)T^{p). (7) 

After projection onto the two Dirac amplitudes we have 



A{p') = l-g'-tr 



p2 



(8) 



B{p') = m,m + Q'\tr [7^ S{p) r,{p)] . (9) 

Equations ([8]) and ([9]) are solved simultaneously with a vertex function calculated at a 
specified order n of vertex dressing. 

The calculations of the Dirac algebra for the vertex diagrams was performed using 



computer-algebraic methods provided by the FeynCalc package ^^;] for Mathematica |16 |. 
The color factors were calculated by numerical contraction of the SU(3) color matrices. The 
number of possible diagrams grows extremely fast with the number of gluon lines included, 
making it difficult to advance too far with the number of loops. For example, while there 
are only 4 possible diagrams with exactly 2 internal gluon lines, there are 27 with 3 and 
38, 232 with 6 internal lines. Such a rapid increase in the number of possible contributions 
forced us to use parallel computing to calculate the invariant amplitude functions for the 
vertex with exactly 6 gluon lines within a reasonable time frame. The contributions to 
include more than 5 x 10^ possible diagrams, which we did not consider feasible to calculate. 
The comparison of the number of diagrams with the previous dressing schemes is shown in 
Table E 

The numerical solutions of the gap equation constructed with the full vertex are shown 
in Figs. [6] and [7] for a current quark mass m = 17 MeV. No attempt has been made to 
simulate 3-gluon coupling - i.e., the color factors used are as directly given by the SU{3)c 
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TABLE I: Comparison of the number of diagrams with exactly n gluon hnes included in various 
vertex dressing schemes: ladder rainbow summed 9|, self-consistently improved ladder rainbow 



summed 



10|, and the present full vertex obtainable with just 2-point gluon lines. 



n 


LR Summed 


Improved 


Full 


2 


1 


3 


4 


3 


1 


(3 


27 


6 


1 


21 


38, 232 


7 


1 


28 


~ 5 * 10^ 



algebra. As found in previous analysis with this model, the results for even and odd n behave 
differently. The even orders provide a behavior for A{s) and M{s) that is analytic and well- 
defined in the infrared spacelike region continuing into the timelike region s > — 1 explored 
here. Mass generation is strong and investigations of the approach to the chiral limit m —>■ 
indicate that these solutions are continuously connected to the Nambu-Goldstone mode of 
spontaneous breaking of chiral symmetry. 

The odd n solutions n > 3 develop non-analytic points in the infrared. There is also a 
tendency for significantly weaker mass generation in such solutions. Evidently, repulsive 
color factors dominate for odd n more than for even n and the mass generation effects are 
in accord with that. 

There is no convincing evidence of convergence with respect to increasing n. At least 
two factors are at work here. Firstly, with the increasing number n of gluon lines, there 
are diagrams at each order that have color factors that are comparable to that of the n = 1 
contribution to the vertex. The class of diagrams omitted in previous studies of this type, 
e.g.. Fig. [T]i, contains members that have larger color factors than the improved ladder- 
rainbow class. Secondly, the number of diagrams specified by n increases rapidly with n, as 
illustrated in Table [B 



2. Effect of 3-point gluon function dressing 

Explicit calculations with gluon 3-point functions proved to be extremely difficult, both 
from the point of view of reproducing all possible diagrams at some non-trivial order in 
quark-gluon coupling and from the fact that no suitable, model gluon 2-point function is 

10 




FIG. 6: (Color Online) Quark propagator amplitude A{s) versus Euclidean s = p^. We use the 
interaction mass scale = 1 GeV and the current mass is m = 0.0175 Q = 17.5 MeV. We show 
the influence of vertex dressing to order n in the number of gluon lines as described in the text, 
n = yields the solid curve and the result is the ladder-rainbow truncation. The other curves are 
n = 1 (long dashed curve, 1-loop vertex), n = 3 (short dashed curve), n = 5 (dot short dashed 
curve) and n = 6 (dot long dashed curve). 

known to us which would allow for feasible calculations beyond one or two loops. Hence we 
resort to the method introduced in Ref. 9|, which allows one to effectively account for 3- 
point gluon function dressing of the vertex through the introduction of the phenomenological 
parameter C, which for iVg = 3 can take values in the range —1/8 < C < 1, as discussed 
after Eq. 

We chose to implement this phenomenological scheme only for the sub-class of vertex 
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FIG. 7: (Color Online) Quark mass function M{s) versus Euclidean s = p^. We use the interaction 
mass scale G = 1 GeV and the current mass is m = 0.0175 G = 17.5 MeV. We show the influence 
of vertex dressing to order n in the number of gluon lines as described in the text, n = yields 
the solid curve and the result is the ladder-rainbow truncation. The other curves are n = 1 (long 
dashed curve, 1-loop vertex), n = 3 (short dashed curve), n = 5 (dot short dashed curve) and 
n = 6 (dot long dashed curve). 

diagrams that correspond to the improved ladder structure where its motivation is unam- 
biguous. We note that all such diagrams have the same color factors as illustrated in Eq. (jS]). 
The results show that for the parameter range C G (0.375, 0.8) the solutions of the gap equa- 
tion have significant mass generation even for the odd n values. This is indicative of the 
attractive effect of the 3-gluon coupling mechanism. Also, in this circumstance there is a 
much weaker tendency of the odd orders > 3 to develop non-analytic points in the deep 
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FIG. 8: (Color Online) Quark propagator amplitude A{s) versus Euclidean s = p^. We use the 
interaction mass scale Q = 1 GeV, the current mass is m = 0.0175 Q = 17.5 MeV. Partial 
account has been taken of 3-gluon coupling through the parameter C = 0.5 which induces more 
mass generation. The influence of quark-gluon vertex dressing to order n in the number of gluon 
lines: n = (solid curve, ladder rainbow truncation), n = 3 (long dashed curve), and n = 5 (short 
dashed curve). The dot dashed curve shows the Pade approximation to the vertex. 

infrared. The solutions for C = 0.5 are presented in Figs. [8] and [91 

3. Pade Approximant 

One can see from Figs. (IE]) and ([9]) that even with the enhanced mass generation generated 
by the strong attractive role of the effective 3-point gluon coupling, the solutions do not show 
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FIG. 9: (Color Online) Quark propagator amplitude M{s) versus Euclidean s = p^. We use the 
interaction mass scale Q = 1 GeV, the current mass is m = 0.0175 Q = 17.5 MeV. Partial 
account has been taken of 3-gluon coupling through the parameter C = 0.5 which induces more 
mass generation. The influence of quark-gluon vertex dressing to order n in the number of gluon 
lines: n = (solid curve, ladder rainbow truncation), n = 3 (long dashed curve), and n = 5 (short 
dashed curve). The dot dashed curve shows the Pade approximation to the vertex. 

convincing convergence with respect to the maximum number, n, of gluon lines allowed 
in the quark-gluon vertex dressing, even in the space-like region. To seek more reliable 
information, we employ a Pade approximant to re-sum the quark self-energy contributions 
with respect to n. While the details of this calculation are presented in Appendix A, the 
graphs corresponding to the solutions of the gap equation with Pade approximant used for 
the quark self-energy contribution are shown in magenta in Figs. ([8]) and (Q. It is clear from 
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the graph that our results are significantly different from those obtained using LR truncation 
and illustrate that convergence is still remote for n ^ 5. 

Unfortunately the present model does not allow consistent solutions in the time-like 
region, where the propagators are needed in order to calculate the meson masses by solving 
the corresponding Bethe-Salpeter equations. 

III. SUMMARY 

We have extended recent studies of the effects of quark-gluon vertex dressing upon the 
solutions of the Dyson- Schwinger equation for the quark propagator. The approach we have 
followed uses a momentum delta function to represent the dominant infrared strength of 
the effective gluon propagator so that the resulting integral equations become algebraic. 
When the guark-gluon vertex is constructed from diagrams that involve only 2-point gluon 
lines, the algebraic form allows the summation of all vertex diagrams, self-consistently with 
solution of the quark propagator. Previous studies along these lines have treated the ladder- 
rainbow class of vertex diagram, and subsequently an improved ladder-rainbow class in 
which all internal vertices are self-consistently dressed. Whereas the earlier improvements 
showed only minor differences from the Ladder-Rainbow approximation, the inclusion of all 
two gluon lines at a given order leads to much larger corrections in the infrared space-like 
as well as the time-like region. 

As far as we are able to take the calculations, there is no evident convergence with 
increasing n. In particular, the behavior with n even or odd is rather different. At least 
two factors are at work here. Firstly, with increasing n, there are diagrams at each order 
that have color factors which are comparable with that of the n = 1 contribution to the 
vertex. The class of diagrams omitted in previous studies of this type, contains members 
that have larger color factors than the improved ladder-rainbow class of diagrams. Secondly, 
the number of diagrams that must be included increases rapidly with n. 

The resulting solutions of the gap equation with the corresponding vertex calculated up 
to a given order, n, in gluon lines show significant deviation from those calculated previously 
using either a bare vertex or a vertex dressed according to a ladder-rainbow sum or its self- 
consistent improvement. The investigation of the infrared behavior of the ghost and gluon 
n-point functions of Refs. [isl indicate that the contributions of ghost loops in dressing 
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the quark-gluon vertex are dominant. It is possible that a consistent treatment involving 
2-, 3-, and 4-point gluon functions as well as ghost contributions might enable a convergent 
summation of diagrams for the vertex. However the present simplified gluon propagator 
permits only zero momentum and this is not compatible with an explicit treatment of ghost 
2-point, and 3- and 4-point gluon functions. In spite of the the numerous interesting applica- 
tions of the Munczek-Nemirovsky model, it is also possible that the advantages of the simple 
gluon propagator, in allowing algebraic simplification of the summation of vertex diagrams 
within the quark DSE , are outweighed by model artifacts. It would be very important to 
explore ways to repeat this calculation with a more realistic gluon propagator. 
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APPENDIX A: PADE APPROXIMANT FOR THE QUARK SELF-ENERGY 

The solutions of the gap equation ([7]) for the quark propagator amplitudes A{s) and B{s) 
depend on the number n of the included gluon lines in the calculations through the Dirac 
projections of the quark self-energy term with appropriately dressed quark-gluon vertex 
function on the RHS of Eqs.® and (Q (the terms under the tr). We re-summed the 
perturbative solutions of the gap equation and obtained a solution at n = oo by employing 
a Pade approximant in n for these self-energy terms. 

In order to construct a Pade approximant for each of these terms we established their 
dependence on n by introducing a fictitious coefficient, uj, to the bare quark-gluon vertex 
coupling g in dressing the vertex function and constructing the projections of the self-energy 
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term in Eqs. ([8]) and ([9]). In each case a Pade approximant of the following form 

1 + + ' ^^^^ 

where A = u"^, was used to match the appropriate self-energy projection term by solving for 
the coefficients and 6j. We considered the self-energy projection terms with the quark- 
gluon vertex dressed to order n = 3 in the number of gluon loops to determine all the 
parameters in ( lAip . 

After determining the coefficients and bi, the parameter uj was set to 1 and the cor- 
responding coupled equations for A{s) and B{s) were solved. For the case C = 0.5, the 
solutions obtained with the approximant self-energy projection terms are finite and contin- 
uous in the space-like region. 
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